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ABSTRACT 

We present Infrared Space Observatory (ISO) observations of several OH, CH and H2O rotational 
lines toward the bright infrared galaxies NGC 253 and NGC 1068. As found in the Galactic clouds in 
Sgr B2 and Orion, the extragalactic far-IR OH lines change from absorption to emission depending 
on the physical conditions and distribution of gas and dust along the line of sight. As a result, most of 
the OH rotational lines that appear in absorption toward NGC 253 are observed in emission toward 
NGC 1068. We show that the far-IR spectrum of OH can be used as a powerful diagnostic to derive 
the physical conditions of extragalactic neutral gas. In particular, we find that a warm (Tfe^l50 K, 
n(H2)<5xl0 4 cm -3 ) component of molecular gas with an OH abundance of ~10~ 7 from the inner 
<15" can qualitatively reproduce the OH lines toward NGC 253. Similar temperatures but higher 
densities (~5xl0 5 cm -3 ) are required to explain the OH emission in NGC 1068. 
Subject headings: galaxies: nuclei — galaxies: ISM — galaxies: individual (NGC253, NGC1068) - 
infrared: galaxies — ISM: molecules 



1. INTRODUCTION 

Galactic nuclei play a key role for our understanding 
of galactic evolution. The use of molecular diagnostics 
to extract their physical conditions helps to understand 
the nature of the prevailing scenarios, e.g., starbursts 
vs. active galactic nuclei (AGN). Space observations 
at mid- and far-IR wavelengths offer an unique oppor- 
tunity to observe the pure rotational lines of light hy- 
drides, the simplest chemical building blocks. In par- 
ticular, H2 observations have revealed the presence of 
a considerable amount of "warm" (excitation tempera- 
tures, T ex , around 100-200 K) neutral gas in the Galactic 
Center (GC; Rodriguez-Fernandez et al. 2001) and also 
in galaxies (Rigopoulou et al. 2002). 

ISO has given the opportunity to observe polar species 
that are sensitive to a broad range of densities and ki- 
netic temperatures. The far-IR spectrum, and thus the 
molecular content, of molecular clouds such as Sgr B2 
in the GC (Goicoechea et al. 2004) is surprisingly sim- 
ilar to that observed in ultra-luminous IR galaxies such 
as Arp 220 (Gonzalez-Alfonso et al. 2004). One of the 
species with the largest abundance is the OH radical. 

Bright IR galaxies containing OH usually show mega- 
maser emission at 18 cm. These lines arise from the 
hyperfme splitting of the OH 2 n 3 / 2 J=3/2 ground rota- 
tional state. Among the possible excitation mechanisms 
of these masers; near- and far-IR pumping, collisions 
and UV radiation, the far-IR excitation seems to be the 
dominant mechanism inverting the population of the hy- 
perfme levels (e.g. Skinner et al. 1997). Furthermore, 
the relative population of the OH levels is governed by 
the far-IR pure rotational lines and detailed treatment 
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of the rotational excitation is required to analyze the OH 
radio lines at 18 cm but also the 2 n 3 / 2 J=5/2 and 2 Hi/2 
J=l/2 lines at 5 and 6 cm (Randell et al. 1995). 

The OH ground-state rotational line at ~119 ptm was 
first detected in absorption toward Sgr B2 (Storey et al. 
1981), while several excited rotational lines have been an- 
alyzed by Goicoechea & Cernicharo (2002). A detailed 
analysis of the far-IR spectrum of OH has been also car- 
ried out in the Orion region (Watson et al. 1985; Vicuso 
et al. 1985; Melnick et al. 1987, 1990; Betz & Boreiko 
1989). The observed lines appear in emission and some 
of them even show more complex P-Cygni profiles. Dif- 
ferent excitation conditions determine the distinct OH 
spectral signatures observed toward Orion and Sgr B2. 

These studies have shown that OH can be a powerful 
prove of the warm neutral gas in molecular clouds. Scaled 
to extragalactic studies, OH can thus provide significant 
information about star forming regions and the heating 
processes in the inner arcsecs of a galaxy, in addition to 
the understanding of the megamaser excitation mecha- 
nism. Moreover, OH is a key molecule in the gas-phase 
chemistry networks. The presence of OH in the warmer 
regions may enhance and even dominate the formation of 
other O-bearing species, so that the determination of its 
abundance is also relevant to understand the prevailing 
chemistry. 

In this letter we present and analyze several OH far-IR 
lines, involving levels up to ~420 K, toward the Starburst 
galaxy NGC 253 and the Seyfert 2 galaxy NGC 1068. We 
compare these extragalactic OH lines with our far-IR ob- 
servations of Sgr B2(M) and Orion IRc2, and run simple 
non-LTE radiative transfer models to qualitatively un- 
derstand the behavior of the detected lines. 

2. OBSERVATIONS AND DATA REDUCTION 

Several pure rotational lines of OH appear in the wave- 
length coverage of the Long- and Short- Wavelength 
Spectrometers (LWS: Clegget al. 1996; SWS: de Graauw 
et al. 1996) on board ISO (Kessler et al. 1996). Brad- 
ford et al. (1999) have previously presented LWS Fabry- 



2 



Goicoechea et al. 




OH 





1.02 
1 

0.9B 



H P ? 



Fig. 1. — ISO/LWS grating observations of NGC 253 and NGC 1068. The ordinate corresponds to the flux and the abscissa is the rest 
wavelength in fim. Atomic fine structure lines are labeled. Zooms around selected wavelengths are also shown in the smaller panels with 
models of the ground— state rotational lines of OH, CH and H2O (see text). The ordinate of these panels is the continuum normalized flux. 



Perot (FP) observations of the OH ground-state line at 
-119 /im toward NGC 253 {v LSR ^ 250 km s^ 1 ), while 
Spinoglio et al. (1999) presented the full LWS grating 
spectrum of NGC 1068 (v LSR ~ 1140 km s" 1 ). In 
this work, we have used the public ISO Data Archive 
(IDA) and analyzed the target dedicated time num- 
bers (TDTs) 24701103, 56901708 (for NGC 253) and 
60500401 and 60501183 (for NGC 1068). The spectral 
resolution of the LWS/grating is A/AA-200. The OH 
2 ni/ 2 ~ 2 n 3 /2 J=5/2-3/2 cross-ladder line at —34 /mi 
was observed with the SWS/grating with a spectral res- 
olution of A/AA-1500 (TDTs 37902123 and 64302927). 
The LWS beam has a -80" diameter while the SWS 
aperture in the —34 /im range is 17"x40". LWS and 
SWS data were analyzed using the ISO Spectral Analysis 
Package 3 (ISAP). Typical routines include: deglitching 
spikes due to cosmic rays, oversampling and averaging 
individual scans, removing baseline polynomials and ex- 
tracting line fluxes by fitting gaussians. More details can 
be found in the LWS handbook (Gry et al. 2003). The 
full ISO/LWS spectra of NGC 253 and NGC 1068 are 
shown in Figure 1 while several OH rotational lines are 
shown in Figure 2. 

3. RESULTS 

The far-IR spectrum of NGC 253 and NGC 1068 
(Fig. 1, large panel) is dominated by the thermal emis- 
sion of dust and by atomic fine structure emission lines 
(Negishi et al. 2001). The signal-to-noise ratio of these 

3 ISAP is a joint development by the LWS and SWS Instruments 
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IPAC, MPE, RAL, and SRON. 



spectra is lower than the ISO/LWS spectrum of Arp 220, 
which is very rich in molecular features (see Gonzalez- 
Alfonso et al. 2004). However, the ground-state rota- 
tional lines of OH, CH and H2O are clearly observed in 
NGC 253 and possibly in NGC 1068 (Fig. 1, small pan- 
els). Figure 2 shows in more detail the OH —34, —53 and 
—79 /im cross-ladder lines, the —84 and —119 /im lines 
of the 2 n 3 / 2 ladder, and the —163 /im line of the 2 Hi/2 
ladder toward NGC 253 and NGC 1068. For comparison 
we also show the OH spectra of archetype star forming 
regions in our Galaxy: Sgr B2 and Orion IRc2. It is 
interesting to note that the behavior of the OH lines in 
NGC 1068 follows that of Orion IRc2, i.e., emission in all 
the observed lines (except in the —53 and —34 /im lines). 
On the other hand, the OH lines in NGC 253 have a sim- 
ilar behavior than in Sgr B2(M), i.e., absorption in the 
2 n 3 / 2 ladder and cross-ladder lines and emission in the 
2 Hi/2 ladder (Goicoechea & Cernicharo 2002). Hence, 
the OH rotational lines show distinct trends in two pro- 
totype galaxies where the nuclei activity is also different. 

Due to its large dipole moment and rotational con- 
stant, radiative pumping of the rotational levels must be 
taken into account as long as collisions to compute the 
OH level population. For the typical temperatures of 
molecular clouds, LTE can only be achieved at very high 
H2 densities, n(H2)>10 10 cm -3 . Hence, the knowledge 
of the dust continuum properties and its distribution re- 
spect to the molecular gas is needed to analyze the OH 
excitation. 

4. DISCUSSION 

In order to qualitatively understand the contribution 
of the physical and geometrical parameters to the OH 
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Fig. 2.— ISO observations of OH rotational lines towards Sgr B2(M) (from Goicoechea & Cernicharo 2002), NGC 253, NGC 1068 and 
Orion IRc2. The ordinate scale corresponds to the continuum normalized flux. The abscissa in all the boxes is the rest wavelength in 
microns. The vertical thick lines indicate the rest wavelengths of the different OH A— doublets. The A-doublets are spectrally resolved in 
the Fabry-Perot observations otherwise grating observations are presented. The rotational energy diagram of OH (neglecting hyperfine 
structure) and the observed transitions are also shown in the middle of the figure. The A-doubling has been exaggerated for clarity. 



spectrum, we have modeled the first 20 rotational levels 
of OH with two radiative transfer codes (LVG and non- 
local; Cernicharo et al. 2000). The hyperfine structure 
of OH was not included in the calculations. Collisional 
cross sections of OH with H 2 from Offer ct al. (1994) 
have been used with an H 2 ortho-to-para ratio of 3. 

It is obviously more difficult to model the molecular 
emission/absorption from a galactic nucleus (an ensem- 
ble of molecular clouds) than to model a single molec- 
ular cloud. Since the observed extragalactic OH lines 
are not velocity resolved and the angular resolution is 
larger than the dimensions of the observed sources, we 
have just considered a simplified model of an IR core + 
OH shell as a first approximation. Therefore, we have 
adopted an spherical geometry consisting of a shell of 
radius R s heii and continuum core of radius R CO re- The 
exact radius determine the resulting OH column density 
and the properties of the dust emission, however, ISO 
observations can not provide that information. Then, 
we have compared the synthetic continuum-normalized 
spectrum with the ISO observations considering that the 
source is composed of several clouds with the same physi- 
cal conditions but do not overlap in the line-of-sight (see 
the discussion by Gonzalez-Alfonso et al. 2004). 



In the case of NGC 253, the mid-IR continuum is con- 
fined to the inner <~10" of the galaxy (Telesco et al. 
1993). We have taken this as a lower limit to R CO re in the 
far-IR. In the other hand, low-J CO observations reveal 
that the molecular gas extends up to ^50" (Harrison et 
al. 1999). However, CO J=7-6 (£„~150 K) observa- 
tions show that the bulk of this emission is confined to 
the inner ~15" (Bradford et al. 2003). In addition, sev- 
eral OH masers at ^18 cm have been detected within 
the central ~10" of the nucleus (Frayer et al. 1998). We 
have assumed that the far-IR OH lines and the IR con- 
tinuum arises from the ~15" region, but we note that a 
more detailed model of NGC 253 could include an ex- 
tended component of low excitation OH gas (that will 
enhance the absorption of the ~119 ^m fundamental 
line; see Gonzalez-Alfonso ct al. 2004 for the case of 
Arp 220). In fact, observations of the extended absorp- 
tion of OH 18 cm lines reveal that OH in the ground- 
state is present as far as ~45" from the nucleus (Turner 
1985). However, here we are mainly interested in the 
general behavior of the higher energy OH lines and since 
the ISO/LWS angular/spectral resolution is poor, we will 
not add more complexity to the calculations in this work. 
We have taken R core /R s ^ e H=0.9, although we also inves- 
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tigate the effect of the dilution of the radiation arising 
from the core. From a gray-body fitting of the LWS 
continuum emission we have adopted a dust opacity law 
T\ = rii 9 [119/A(^m)] 1 - 5 with r U9 =l and T d =41 K. 

To compare with the NGC 1068 observations, we have 
considered that the OH emission also arises from a single 
component with similar geometry and continuum emis- 
sion as for NGC 253. However, our calculations showed 
that the OH emission fluxes could only be reproduced 
if the emission arises from a more compact component 
with 9 s heii^0" and 9 core ^5". This suggests that the 
OH emission may be dominated by dense gas in cloud 
cores near the nucleus. This is the case of HCN observed 
at higher resolution (Jackson et al. 1993; Tacconi et al. 
1994). Note that OH masers at —18 cm have been found 
in the inner arcsec of NGC 1068 (Gallimore et al. 1996). 
Higher angular resolution is needed to estimate the con- 
tribution of the Starburst ring to the OH emission. 
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Fig. 3. — Large Velocity Gradient models of different kinetic 
temperature and density for the observed OH lines in NGC 253. 
The OH column density is 6.5 xlO 16 cm -2 . The thick contour 
corresponds to the equivalent temperature of the continuum source 
(Rco7'e/R. s heii = 0-9 and T119 (dust)=l). To the left of this contour, 
lines are predicted in absorption. 

Several results from the LVG computations are 
shown in Fig. 3. The excitation temperature of each 
rotational transition as a function of gas temperature 
(Tfc) and density is shown in each panel. The thick 
black contour represents the equivalent temperature 
of the continuum. Therefore, to the left of this con- 
tour, the OH lines are predicted in absorption. From 
the excitation models it is clear that high densities 
are required to observe the OH lines in emission, 
i.e., the case of NGC 1068. If the R C ore/Rsheii ratio 
is decreased or the shell thickness is increased, the 
iso-T ea; contours on Fig. 3 (the absorption/emission 
boundary), will move to the left toward lower densities. 
This is due to the increase of the reemission from 
the gas that do not absorb the core emission. Note 
that reemission also scales with the thickness of the shell. 

To investigate the effect of the radiative coupling be- 
tween components of different physical conditions we 
have also run several nonlocal models (see Gonzalez- 
Alfonso & Cernicharo 1993) for OH including both the 
gas and dust transfer (Cernicharo et al. 2000). We run 
calculations for different geometrical and physical condi- 
tions. As an example, Fig. 4 shows some of the results 



for the -163.12 and -79.18 /xm lines. The 2 n 1/2 J=3/2- 
1/2 line at —163 /jm is observed in emission in all the 
sources. Goicoechea & Cernicharo (2002) showed that 
the —163 /im emission in Sgr B2(M), were the far-IR con- 
tinuum is optically thick, is produced by a fluorescent- 
like mechanism in low density (<10 4 cm -3 ) warm gas 
(T fe ~150-300 K). In the other hand, Melnick et al. (1987) 
showed that T fc <100 K and n(H 2 )>10 6-7 cm -3 are re- 
quired to fit the observed emission in Orion IRc2 shocked 
region, so that collisions dominate the excitation. How- 
ever, these authors argued that even with those condi- 
tions, the radiative pumping to the levels that give rise 
to the —163 /jm emission has to be taken into account. 





-100 100 -100 100 



v rest [ km s '] v rest [ km s '] 

Fig. 4.— Nonlocal models for the 2 n 1/2 - 2 n 3/2 J=l/2+-3/2 _ 
transition at 79.18 (im (upper panel) and for the 2 n 1 y 2 J=3/2+- 
2/2~ transition at 163.12 /im (lower panel). Line intensities have 
been normalized to the continuum emission. The abscissa is the 
rest velocity. The velocity resolution in the line profiles is 1 km s — 1 . 
All the models have x(OH)=10 — 7 . Left panels: Models with 
R CO r e /R s heii=0.9, T fc =150 K, n(H 2 )=5xl0 4 cm -3 and variable 
dust opacity of the continuum source. Right panels: Models with 
Rcore/Rsheii = 0-9i Tng(dust)=l and variable Tj, and n(H 2 ). 

From our calculations it is clear that while at low 
OH abundances [x(OH)] the —79, —53 and —34 /zm 
cross-ladder absorption depths are proportional to the 
OH column density [A(OH)], the —163 /im emission 
over the continuum flux is nearly constant because the 
emission mechanism at low densities depends mainly on 
the continuum radiation. At high x(OH), reemission 
start to dominate the —79 and —119 /im lines. Taking 
into account the simplicity of our calculations, we find 
x(OH)— 10~ 7 as the most realistic value for NGC 253. 
Lower values are required if the core emission is more 
diluted, otherwise reemission dominates. On the other 
hand, larger abundances can be possible if OH emis- 
sion/absorption and far-IR continuum emission have dif- 
ferent dilution factors. 

Density and T& determine the role of collisional excita- 
tion. For low n(H2), radiative excitation dominates, and 
the OH spectrum is less sensitive to n(H.2) and T& vari- 
ations. This seems to be the case of NGC 253 inner re- 
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gions and the GC molecular clouds, which probably share 
similar physical conditions. The effect of Tfe is more rel- 
evant for higher dilution of the core radiation. We found 
that the model with T fe -150 K and n(H 2 )<5xl0 4 cm" 3 
better fits the OH observations of NGC 253. This tem- 
perature is in between the T ex derived from the H 2 5(0) 
line by Rigopoulou et al. (2002) and the T& derived from 
the CO J=7-6 line emission by Bradford ct al. (2003). 
Finally, from the ground-state rotational lines of CH and 
H2O toward NGC 253 and assuming intrinsic line widths 
of ^200 km s _1 and T ex <Td (absorption lines) we found 
a lower limit for the CH and H 2 column densities of 
~5xl0 13 and ~lxl0 15 cm~ 2 respectively (Fig. 1). 

Similar arguments also apply for NGC 1068. Fol- 
lowing Sternberg et al. (1994) we have first assumed 
Tfc=50 K. Taking into account the simplifications in 
our models, the calculations for NGC 1068 show that 
the density must be larger than 10 6 cm~ 3 , otherwise 
the OH lines are still predicted in absorption. For 
ra(H 2 )=5xl0 5 cm~ 3 , the temperature has to be in- 
creased to ~150 K to fit the observations, so that it is 
likely that a fraction of warm neutral gas also exists 
in the clouds near the nucleus of NGC 1068. Intense 
emission in other OH rotational lines (not detected by 
ISO) are predicted for higher densities and temperatures. 
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5. CONCLUSIONS 

Irrespective of the specific nuclei scenario and of the 
density of the dominant molecular clouds, OH obser- 
vations show that a similar component of warm neu- 
tral gas (T fc -150 K) exists in NGC 253 and NGC 1068. 
This component can influence the chemistry as neutral- 
neutral reactions involving OH start to be efficient for 
the production of other oxygen molecules such as H 2 
and SiO. The contribution of PDRs, shocks, X-rays and 
cosmic rays to the heating of such component in Seyfert 
vs. Starburst galaxies is still far from settled. Future 
works using diagnostics of the warm gas will be extremely 
useful. In particular, several rotational lines of light hy- 
drides such as OH, CH and H 2 will be observed at 
larger angular resolution by the HIFI and PACS instru- 
ments, on board Herschel in many galactic and extra- 
galactic sources. 
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